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Abstract 



o 

. We prove general nonlinear stability and existence theorems for rotating star solu- 

tions which are axi-symmetric steady-state solutions of the compressible isentropic Euler- 
Poisson equations in 3 spatial dimensions. We apply our results to rotating and non- 
rotating white dwarf, and rotating high density supermassive (extreme relativistic) stars, 
^ • stars which are in convective equilibrium and have uniform chemical composition. This 

bX), paper is a continuation of our earlier work f[28|). 
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1 Introduction 

The motion of a compressible isentropic perfect fluid with self-gravitation is modeled by the 
Euler-Poisson equations in three space dimensions (cf [6]): 

pt + V- (pv) = 0, 

(pv)t + V • (pv v) + Vp(/)) = -/9V$, (1.1) 
A<1) = Airp. 
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Here p, v = {vi,V2,V3), p{p) and ^> denote the density, velocity, pressure and gravitational 
potential, respectively. The gravitational potential is given by 



where * denotes convolution. System is used to model the evolution of a Newtonian 
gaseous star ([6j). In the study of time- independent solutions of system (jl.ip . there are two 
cases, non-rotating stars and rotating stars. An important question concerns the stability of 
such solutions. Physicists call such star solutions stable provided that they are minima of an 
associated energy functional ([37], p. 305 & [33]). Mathematicians, on the other hand, consider 
dynamical nonlinear stability via solutions of the Cauchy problem. The main purpose of this 
paper is to prove a general theorem which relates these two notions and shows that for a wide 
class of Newtonian rotating stars, minima of the energy functional are in fact, dynamically 
stable. This is done for various equations of state p = p{p) which includes polytropes, 
supermassive, and white dwarf stars. 

For non-rotating stars, Rein ([32]) has proved nonlinear stability under various hypotheses 
on the equation of state, including in particular, polytropes where p = kp"^, 7 > 4/3; his theory 
applies to neither white dwarf nor supermassive stars. In a recent paper, [28], we studied 
nonlinear stability of rotating polytropic stars, where p = kp'^ , 7 > 4/3. In this paper, we 
generalize these results to rotating white dwarf and supermassive stars, thereby completing 
the nonlinear stability theory for rotating (and non-rotating) compressible Newtonian stars*. 

Our main theorem applies to minimizers of an energy functional with a total mass con- 
straint. The crucial hypotheses are that the infimum of the energy functional in the requisite 
class, be finite and negative. This is verified for both white dwarf and supermassive stars by 
combining a scaling technique used by Rein ([31]), together with our method in [28] where we 
use some particular solutions of the Euler-Poisson equations in order to simplify the energy 
functional. It should be noticed that neither the scaling technique in [3T] nor the method 
in |28j using particular solutions of Euler-Poisson equations apply to white dwarf stars di- 
rectly. As a bi-product of our method, we prove the existence of a minimizer for the energy 
functional, which is a rotating white dwarf star solution, in a class of functions having less 
symmetry than those solutions obtained in [2] and [13] . The method in [2] and [13] is to con- 
struct a specific minimizing sequence of the energy functional, each element in the sequence 
being a steady solution of the Euler-Poisson equations. In contrast, our method is to show 
that any minimizing sequence of the energy functional must be compact (cf. Theorem 3.1 
below). This fact is crucial for both existence and stability results. 



* In all cases under consideration, stability is only "conditional" because no global in time solutions have 
been constructed so far for compressible Euler-type equations in three spatial dimensions; this is a major open 
problem. 




(1.2) 
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For a white dwarf star (a star in which gravity is balanced by electron degeneracy pressure) , 
the pressure function p{p) obeys the following asymptotics ([6], Chapter 10): 

P{p) = Cl//^ - C2P^/^ + • • • , CO, 
p(p) =dlp5/3_d2pV3 + (9(^3)^ ^^0, 

where ci, C2, di and d2 are positive constants. The existence theory for non-rotating white 
dwarf stars is classical provided the mass M of the star is not greater than a critical mass 
Mc {M < Mc) ([6j). For rotating white dwarf stars with prescribed total mass and angular 
momentum distribution, Auchumuty and Beals ([2]) proved that if the angular momentum 
distribution is nonnegative, then existence holds if M < Mc- Friedman and Turkington 
([H]) proved existence for any mass provided that the angular momentum distribution is 
everywhere positive; see Li ([23j), Chanillo & Li ([7]) and Luo & Smoller (23) for related 
results for rotating star solutions with prescribed constant angular velocity. To the best of 
our knowledge, our stability theorem in this paper for rotating and non-rotating white dwarf 
stars with M < M^. is the first nonlinear dynamical stability theorem for such stars. 

For a supermassive star (a star which is supported by the pressure of radiation rather than 
that of matter; sometimes called an extreme relativistic degenerate star [33]), the pressure 
p{p) is given by ([37]): 

p{p) = kp\ 7 = 4/3, (L4) 

where /c > is a constant. For non-rotating spherically symmetric solutions for supermassive 
stars, Weinberg ([37j) showed that the total energy vanishes; thus to quote Weinberg ([37]. 
p. 327) "the polytrope with 7 = 4/3 is trembling between stability and instability", and 
he remarks that one needs to use general relativity to settle this stability problem. For 
rotating supermassive star solutions, we show here that the energy is negative £^ < due to 
the rotational kinetic energy (see (4.26) below). Thus the stability problem falls within the 
framework of Newtonian mechanics and so our general stability theorem applies to show that 
rotating supermassive stars are nonlinearly stable, provided that M < Mc- 

For the stability of both white dwarfs and supermassive stars, we require that the total 
mass of each one lies below to a corresponding critical mass, a "Chandrasekhar" limit. We 
show that this holds because the pressure function for both is of the order p^^^ as p — > cxd. 

The above dynamical stability results for rotating stars apply for axi-symmetric perturba- 
tions. For non-rotating stars, G. Rein ([32]) proved nonlinear dynamical stability for general 
perturbations. However, his result does not apply to white dwarf stars. For non-rotating 
white dwarf stars, the problem was formulated by Chandrasekhar [5] in 1931 (and also in 
|llj and [19]) and leads to an equation for the density which was called the " Chandrasekhar 
equation " by Lieb and Yau in [21] . This equation predicts the gravitational collapse at some 
critical mass ([5j and [6j). This gravitational collapse was also verified by Lieb and Yau ( (24j ) 
as the limit of Quantum Mechanics. In Section 5, we prove the nonlinear dynamical stability 
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for non-rotating white dwarf stars with general perturbations provided that the total mass is 
below some critical mass. 

Other related results besides those mentioned above for compressible fluid rotating stars 
can be found in [3], 0], [l2], and [ST]. 

The linearized stability and instability for non-rotating and rotating stars were discussed 
by Lin ([25] ), Lebovitz ([21]) and Lebovitz & Lifschitz ([22]). Related nonlinear stability 
and instability results for gaseous stellar objects can be found in Guo k. Rein ([l5], [16]) and 
Jang ([H]). Related results for the Euler-Poisson equations of self-gravitating fluids can be 
found in [TO], [IS], [22] and 



2 Rotating Star Solutions 

We now introduce some notation which will be used throughout this paper. We use j to 
denote /jjs, and use || • \\q to denote || • ||l<7(ir3). For any point x = {xi^X2^x^) G M^, let 



r{x) = + xl, z{x) = X3, Br{x) ={yG M^ \y - x\ < R}. (2.1) 
For any function / G L^(M^), we define the operator B by 

Bf{x)= [ JM^dy = f*-^. (2.2) 

J \x-y\ \x\ 

Also, we use V to denote the spatial gradient, i.e., V = = {dx^, dx^, dx^)- C will denote 
a generic positive constant. 

A rotating star solution (p, v, $)(r, z), where r = y^xf + X2 and z = X3, x = (xi, X2, 2:3) G 
M^, is an axi- symmetric time-independent solution of system (jl.ip . which models a star ro- 
tating about the xs-axis. Suppose the angular momentum (per unit mass), J{mp{r)) is 
prescribed, where 

rripir) = l p{x)dx = j 27rs / p{s, z)dsdz, (2.3) 



V^xf+x|<r Jo 



is the mass in the cylinder {x = (xi,X2,X3) : \J x\^ x\ < r}, and J is a given function. In 
this case, the velocity field v(x) = (^1,^2,^3) takes the form 

^,(^\ _ ( X2J{mp{r)) xiJ{mp{r)) 

Substituting this in (jl.ip . we find that p{r,z) satisfies the following two equations: 

drp{p) = pdr{Bp) + pL{mp{r)r-^, 
dzP{p) = pdziBp), 
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where the operator B is defined in ()2.2p . and 

L{mp) = J^{mp) 
is the square of the angular momentum. We define 

A{p) = p r^-^ds. (2.5) 
It is easy to verify that (cf. [2]) (|2.4p is equivalent to 

A {p{x)) + / L{mp{s)s ■^ds - Bp{x) = A, where p{x) > 0, (2.6) 

Jr{x) 

for some constant A. Here r{x) and z(x) are as in (j2.ip . Let M be a positive constant and 
let Wm be the set of functions p defined by , 



Wm ={p : ^ M, p is axisymmetric, p > 0, a.e., 

' p{x)dx = M, [ ( A{p{x)) + ^^''^^^^^'1^^''^^ + p{x)Bp{x) ] dx < +00.} 



For p G Wm, we define the energy functional F by 



1 p{x)L{mp{r{x))) 1 
r(x)^ 



F(p) = / [A{p{x)) + - (2.7) 



In (j2.7p . the first term denotes the potential energy, the middle term denotes the rotational 
kinetic energy and the third term is the gravitational energy. 

For a white dwarf star, the pressure function p{p) satisfies the following conditions: 



lim 



0, lim —J- = p{p) > as /9 > 0, 



/3^0+ p4/3 ' p^oo p4 

where is a finite positive constant. Assuming that the function L G C^[0, M] and satisfies 

L(0) = 0, L(m) > 0, for <m< M, (2.9) 

Auchmuty and Beals (cf. [2]) proved the existence of a minimizer of the functional F{p) in 
the class of functions Wm,s = Wm H Ws, where 

Ws = {p:R^^ M, /j(xi, xa, -X3) = p{xi,X2,X3), x,eR,i = l, 2, 3}. (2.10) 

Their result is given in the following theorem. 

Theorem 2.1. (12]). If the pressure function p satisfies 112. 8\) (for either < < +00 or 
K = +00 and \2. 9) holds, then there exists a constant Mc > depending on the constant ^ 
in \2.8\) fif ^= +00 then Mc = +00, ifO<K< +00, then < Mc < +00) such that, if 

M < Mc, (2.11) 
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then there exists a function p{x) G Wm,s which minimizes F{p) in Moreover, if 

G = {xeM.^ : p{x) > 0}, (2.12) 

then G is a compact set in , and p G C^{G) n C^(M^) for some < /3 < 1. Furthermore, 
there exists a constant p, < such that 

A'{p{x)) + J^^^ L{mp{s)s-^ds - Bp{x) = p, xeG, 
l^^^ L{mii{s)s-^ds - Bp{x) >p, xeM.^-G. 

Remark 1. When < < cxd, the constant < Mc < +oo in (j2.1ip is called critical mass. 
The critical mass was first found by Chandrasekhar (cf . [S] ) in the study of non-rotating white 
dwarf stars. When < < oo, it was proved by Friedman and Turkington (|13|) that, if the 
angular momentum satisfies the following condition 

J £ C^([0,M]), J'(m) > 0, for < m < M, J(0) = 0, J(m) > for < m < M, (2.14) 

where J is the angular momentum, then the condition (j2.11|) can be removed, i.e., the above 
theorem holds for any positive total mass M. 

In this paper, we are interested in the minimizer of functional F in the larger class Wm- 
By the same argument as in [2] , it is easy to prove the following theorem on the regularity of 
the minimizer. 

Theorem 2.2. Suppose that the pressure function p satisfies: 

lim 4?^ = 0, lim 4^ = oo, p{p) > as p > 0, (2.15) 

and the angular momentum satisfies i2.9\) . Let p be a minimizer of the energy functional F 
in Wm and let 

r = {x e : p{x) > 0}, (2.16) 
then p £ G{W^) H C^(r). Moreover, there exists a constant A such that 

^'(Pi^)) + IZ) L{m~,{s)s-^ds - Bpix) = A, xGT, 
l^^^ L{mp{s)s-^ds - Bp{x) > A, x G _ p. 

We call such a minimizer p a rotating star solution with total mass M and angular momentum 

3 General Existence and Stability Theorems 

For the angular momentum, besides the condition (j2.9p . we also assume that it satisfies the 
following conditions: 

L{am) > a'^/'^Lim), 0<a<l, 0<m<M, (3.1) 
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L'{m) > 0, < m < M. (3.2) 
Condition (|3.2p is called the Solberg stability criterion ([35j). 

3.1 Compactness of Minimizing Sequence 

In this section, we first establish a compactness result for the minimizing sequences of the 
functional F . This compactness result is crucial for the existence and stability analyses. 

Theorem 3.1. Suppose that the square of the angular momentum L satisfies 112. 9\) . i3. 1\) and 
13. and the pressure function p satisfies the following conditions 

peC^[0,+oo), [ ^dp<+oo,lim^ = K, p{p)>0,p'{p)>Ofoi p>0, (3.3) 
Jo P P-""^ P'^ 

where < K < +oo and 7 > 4/3. If 

(1) 

inf F(p)<0, (3.4) 

P&Wm 

and 

(2)forp£WM, 

I ^ l p(.)L(^ (Kx))) ^ ^ 
J 1 r\xY 

for some positive constants Ci and C2, then the following hold: 

(a) If {p^} C Wm is a minimizing sequence for the functional F, then there exist a sequence 
of vertical shifts OiBs (ai G M, es = (0,0,1)^, a subsequence of {p^}, (still labeled {p^}), and 
a function p G Wm, such that for any e > there exists R > with 

/ Tp\x)dx < e, i G N, (3.6) 
J\x\>R 

and 

Tp'^{x) p, weakly in L'^{M.^), as i ^ 00, (3.7) 

where Tp'-{x) := p^{x + fljea). 

Moreover 

(b) 

VB{Tp') VB{p) strongly in L^iR-^), as i ^ 00. (3.8) 
(c) p is a minimizer of F in Wm ■ 

Thus yO is a rotating star solution with total mass M and angular momentum \/L. 
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Remark 2. i) The assumption ()3.4p is crucial for our compactness and stability analysis. The 
physical meaning of this is that the gravitational energy, the negative part of the energy 
F, should be greater than the positive part, which means the gravitation should be strong 
enough to hold the star together. In section 4, we will verify this assumption. Roughly 
speaking, in addition to (|3.3p . if we require 

lim ^ = a, (3.9) 

for some constants 71 > 4/3 and < a < +00, then (|3.4p holds for the following cases: 

(a) when 7 = 4/3 (where 7 is the constant in (j3.3p ). if the total mass M is less than a "critical 
mass" Mc, then (j3.4p holds. This case includes white dwarf stars. For a white dwarf star, 
71 = 5/3. 

(b) When 7 > 4/3, (|3.4p holds for arbitrary positive total mass M. This generalizes our 
previous result in [28j for the polytropic stars with p{p) = , j3 > 4/3. 

It should be noted that (j3.9p does not apply to suppermassive star, i.e. p{p) = kp^/^ . For 
the supermassive star, in order that have that (j3.4p hold, in additional to requiring that the 
total mass is less than a "critical mass", we also require that the angular momentum (per 
unit mass) J is not identically zero. 

ii) Assumption (2) in the above theorem implies that the functional F is bounded below, i.e., 

inf F(p)>-oo. (3.10) 

We will verify this assumption in Section 4 (see Theorem 4.1). 

iii) The inequality (j3.6p is crucial for the compactness result (j3.8p . One of the difficulties in 
the analysis is the loss of compactness because we consider the problem in an unbounded 
space, M^. The inequality (j3.6p means the masses of the elements in the minimizing sequence 
Tp^{x) "almost" concentrate in a ball Bji{0). 

iv) It is easy to verify that the functional F is invariant under any vertical shift, i.e., if 
p(-) G Wm, then p{x) =: p{x + ae^) G Wm and F{p) = F{p) for any a G M. Therefore, if 
{/?*} is a minimizing sequence of F in Wm, then {Tp*} =:= p'^{x + 0^63) is also a minimizing 
sequence in Wm- 

Theorem 13.11 is proved in a sequence of lemmas with some modifications of the arguments 
in |28j . We only sketch the proofs of those lemmas and Theorem 13.11 Complete details can 
be followed as in [2H]. We first give some inequalities which will be used later. We begin with 
Young's inequality (see [17], p. 146.) 

Lemma 3.1. // / G n U' , l<p<q<r< +00, then 

,<ll/lPl/ll^^ ^ = l~-i~S^ - (3-11) 
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The following two lemmas are proved in [2]. 
Lemma 3.2. Suppose the function f E L^iR-^) n L'?(M3). If 1 < q < 3/2, then Bf =: f 



is in UiW) for3<r < 3g/(3 - 2q), and 

||i?/||.<c(||/||?||/||^-'' + ||/||5||/||J-^), (3.12) 

for some constants C > 0, < b < 1, and < c < 1. If q > 3/2, then Bf{x) is a hounded 
continuous function, and satisfies 1^3. with r = oo. 

Lemma 3.3. For any function f G L^{R^) D L'^/^{M.^), if j > 4/3, then VBf G L'^{M.^). 
Moreover, 

1 1 f{x)Bf{x)dx\ = ^\\VBf\\l<c(^j \f\"\x)dx^ (^j \f\{x)dxY\ (3.13) 

for some constant C . 

We also need the following lemma. 

Lemma 3.4. Suppose that the pressure function p satisfies 113. 3\) and 113. 5\) holds. Let {/)*} C 
Wm be a minimizing sequence for the functional F. Then there exists a constant C > such 
that 

[ [(pT(x) + ^ ^^(^)^Kf(^))) ],, < c, for all . > 1, (3.14) 
J 2 r{xy 

where 7 > 4/3 is the constant in i3.3\j . So, the sequence is bounded in L'''(IR^). 
Proof. By ()3.5p . we know that 

f ^ l p\x)L{rn^.{r{x))) ^ . ^ ^ 

J 2 r{xy 

for any minimizing sequence {p*} C Wm for the functional F, where we have used that 
{F(/3*)} is bounded from above since it converges to inf^y^^^ F. It is easy to verify that, by 
virtue of ([331) and ([23]) . 

A(o) K 

lim — ^ = -, A{p) > for p > 0. (3.16) 

p^oo pi 7 — 1 

Therefore, there exits a constant p* > such that 

aA{p)>p^, forp>p*, (3.17) 
where a = Hence, for p G Wm-, 

p^dx< I {p*y^pdx + a f A{p)dx 
J p<p* J p>p* 

< {p*y-^M + a j A{p)dx. (3.18) 



Applying this inequality to p*, we conclude that the sequence {p*} is bounded in L'^(M ) by 
using (|3T5]) . □ 
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For any M > 0, we let 



fM= inf Fip). (3.19) 

P&Wm 



Lemma 3.5. // holds, then fj^ > {M/Mf/^fM for every M > M >0 . 

Proof. The proof follows from a scaling argument as in [31J and [28j. Take a = [M/Myl^ 
and let p{x) = p{hx) for any p £ Wm- It is easy to verify that p £ Wf^. Moreover, for r > 0, 
it is easy to verify, (as in [S^) that 

rripir) = -^nipiar). (3.20) 
Since L satisfies p.ip and a > 1, we have 



L{mp{r)) > ^L{m,{br)). (3.21) 



Thus, as in |28], we can show that 



p{x)L{mp{r{x))) ^ ]_ f p{x)L{mp{r{x))) 
r(x)^ J r{xY 



Therefore, since a > 1, it follows from (I3.2ip and ()3.22h that 
Fip) > a-3 / Aip)dx - ^ / pBpdx + ^ / 



-5f f ^^ f T. , 1 f p{x)L{mp{r{x))) 

a'^J A{p)dx - - y pBpdx + -J ^ 

{M/Mf/^F{p). (3.23) 



Since p p \s one-to-one between Wm and W^, this proves the lemma. 



□ 



Lemma 3.6. Let C Wm be a minimizing sequence for F. Then there exist constants 
rg > 0, 5o > 0, zq £ N and x' G M'^ with r(x*) < tq, such that 

p\x)dx >6q, i> io. (3.24) 

Proof. First, since limi^oo F{p^) — > /m and /m < (see (3.4)), for large i, 

- ^ < -F{p') <\j p'Bp^dx. (3.25) 

For any i, let 



5i = sup / p\y)dy. (3.26) 

a;eR3 J|y-x|<l 
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Now 



p'Bp\x)dx (3.27) 



[ P\x){ [ + [ + [ }- 

JR^ J\y-x\<l Jl<\y-x\<r J\y-x\>r 



-dydx 



-x\>r \y ^\ 

=:Di + D2 + D3, (3.28) 

and D3 < M^r^^. The shell 1 < |y — x| < r can be covered by at most Cr^ balls of radius 1, 
so D2 < CMdir'^. By using Holder's inequality and applying (|3.12p to the restriction of to 
{y ■.\y — x\ < 1}, we get 



Di < \\p%/3\\ [ 
Jh 



p\y) 



\y-x\<i \y ^1 

^ ri\\ *ll / II i||l— b i ii *I|C|| inl— c\ 

< (-^\\P II4/3 \ \\Xb^{x)P IIiIIp II4/3 + \\Xb^(x)P II1II/5 II4/3 j 

< Wh/, [5\\\p%-j^ + 5t\\p%-j^) , (3.29) 

where < & < 1 and < c < 1. Now since {||/5*||7} is bounded, it follows that {H/O^IU/s} is 
bounded due to the fact 7 > 4/3 in view of ([3TT]) and ||/9*||i = M; this gives Di < C{8\ + 51). 
It follows that we could choose r so large that the above estimates give / p'^Bp^{x)dx < — /m 
if5i were small enough. This would contradict (|3.25|) . So there exists 5o > such that 6i > 6q 
for large i. Thus, as i is large, there exists x* G and £ such that 

p\x)dx > Jo, i> io- (3.30) 

Biix^) 

We now prove that there exists rg > independent of i such that x^ must satisfy r(x*) < tq 
for i large. Namely, since p^ has mass at least 5q in the unit ball centered at x*, and is axially 
symmetric, it has mass > Cr(x*)5o in the torus obtained by revolving this ball around xs-axis 
(or z-axis). Therefore r(x*) < {CSq)~^ M. □ 

In order to prove Theorem 13. H we will need the following lemma. 

Lemma 3.7. Let {/*} be a bounded sequence in L'~'{W^) (7 > 4/3j and suppose 

f f weakly in L^(M^). 

Then 

(a) For any R > 0, 

yB{xBn{^)h ^ VS(xb^(0)/°) strongly in L^M.^), 
where x the indicator function. 

(b) If in addition {/*} is bounded in L^{R^), f^ G ^-'^(M^), and for any e > there exist 
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R > and G N such that 

\f{x)\dx < e, i> io, (3.31) 



'\x\>R 
then 

VBf VBf strongly in L'^{R^). 

Proof. This lemma follows easily from the proof of Lemma 3.7 in |31j . due to the following 
observation: 

The map: p € L'>'(M^) ^ /b^(o)VB(/s^(o)/o) is compact for any i? > 0, if 7 > 4/3. □ 

With above lemmas, the proof of Theorem 13. II is similar to that in [28j. So we only outline 
the main steps. 

Proof of Theorem \3.1\ 

Step 1. We begin with a splitting as in [31j. For p £ Wm, for any < i?i < R2, we have 

P = PX\x\<Ri + PXRi<\x\<R2 + PX\x\>R2 ='■ Pi+ P2 + P3, (3.32) 



where x is the indicator function. It is easy to verify that 

f p{x)L{mp{r{x)) = V / Pj(x)L(mp^. (r(x)) 



r'^(x) 

^sp f Pj(3;)(L(mp(r(x)) - L{mp^{r{x)) 



3 



/ Pi(x)L(mo.(r(x)) 
In the last inequality above, we have used (3.2). So, we have 

3 

F{p)>J2PiPj)- E (3.34) 

j=l l<i<j<3 

where 

If we choose R2 > 2Ri in the splitting ()3.32p . then 

/i3 < (3.35) 

By (;H:T^ and (fmi) . we have 

-^12 + -^23 

= ^ y V(5pi + Spa) • VBp2dx < C\\V{Bpi + Bp3)\\2\\VBp2 

||J/^||Vi?p2||2 < CM'/^p\\lf^\\VBp2\\2. (3.36) 

12 
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Using Lemma [331 p.4p . (j3.34p . (j3.35p and (j3.36p . and following an argument as in the proof 
of Theorem 3.1 in [31], we can show that 



/m - F{p) 

< (1 - {^f" - {^f" - i^fhfM + C{R,^ + M^^M^WVEp^W,) 

< C/mM^Ms + C(i?2 ' + Mi/3||p||2/3||vSp2||2), (3.37) 

by choosing i?2 > 2Ri in the splitting (|3.32p . where Mi = f pi{x)dx {i = 1,2,3.) Let {p*} be 
a minimizing sequence of F in Wm- By Lemma |3.6[ we know that there exists io G N and 
6o > independent of i such that 

/ p\x)dx > 5o, if i > io, (3.38) 

where Oj = z{x^) and Rq = rQ + 1, x* and ro are those quantities in Lemma [3?6l es = (0, 0, 1). 
Having proved (|3.38p . we can follow the argument in the proof of Theorem 3.1 in [31j to verify 
([33T]) for 

f{x) = Tp\x) =:p'(. + a,e3) 

by using (|3.34p and (j3.38p and choosing suitable Ri and R2 in the splitting (I3.32P . We sketch 
this as follows. The sequence Tp^ =: p*(- + a^es), i > io, is a minimizing sequence of F in 
Wm (see Remark 2 after Theorem 13. ip . We rewrite ()3.38p as 

/ Tp\x)dx > (5o, i > io- (3.39) 

Applying (j3.37p with Tp^ replacing p, and noticing that {Tp*} is bounded in L'''(M^) (see 
Lemma 3.4) (so {HT'/O^IU/a} is bounded if 7 > 4/3 in view of (|3.1ip and the fact = M), 

we obtain, if R2 > 2Ri, 

- CfMMiMi < C{R^^ + \\VBTfi\\2) + F{Tp') - /m, (3.40) 

where M{ = jTp\{x)dx = J^^^^^Jp'{x)dx„ Ml = jTp\{x)dx = J^^^^j^^Tp'{x)dx and 
Fp2 = XRi<\x\<R2Fp^- Since {Tp"^} is bounded in L'''(]R^), there exists a subsequence, still 
labeled by {Tp*}, and a function p £ Wm such that 

Tp' p weakly in r'{R^). 

This proves ([321) • By (I339D . we know that M| in (^3lO\\ satisfies M| > 5o for i > io by 
choosing Ri > Rq where Rq is the constant in (j3.39p . Therefore, by (j3.40p and the fact that 
fM <0 (cf. (fOp i , we have 

- CfM^oMl < CR^^ + C\\VBp2\\2 + C\\VBTp^^ - VSpsIb) + F{Tp') - Jm, (3.41) 
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where p2 = X\x\>R2P- Given any e > 0, by the same argument as [3T], we can increase 
Ri > Rq such that the second term on the right hand side of ()3.4ip is smah, say less than 
e/4. Next choose R2 > 2Ri such that the first term is small. Now that i?i and R2 are fixed, 
the third term on the right hand side of (j3.4ip converges to zero by Lemma IS-Tf a). Since 
{T/o*} is a minimizing sequence of F in Wm, we can make F{Tp^) — Jm small by taking i 
large. Therefore, for i sufficiently large, we can make 



Ml 



: [ Tp\x)dx < e. (3.42) 

J\x\>R2 



This verifies (j3.3ip in Lemma 3.7 for /* = Tp'. By weak convergence we have that for any 
e > there exists R > such that 



M -e< p{x)dx < M, 

Jbr(O) 



which implies p G L^(M^) with f pdx = M. Therefore, by Lemma 13.7( b) .we have 

\\SJBTp' -VBp\\2-^Q, i^+00. (3.43) 

This proves dM]). in Theorem O follows from (102]) by taking R = R2. 

Let {p*} be a minimizing sequence of the energy functional F, and let p he a weak limit 
of {Tp*} in L'^(M^). We will prove that p is a minimizer of F in Wm] that is 

F{p) < lim inf F{Tp'). (3.44) 

By ()3.3p . there exist positive constants C and p* such that 

A'{p)<Cp'^~\for p>p*, (3.45) 

where 7 > 4/3 is the constant in (3.3). Since p L"' and j pdx = M, we can conclude 

A'{p) £ L"'' , where L'^' is the dual space of L"', i.e., 7' = In view of (2.5) and (3.3), we 
have 

A"{p) = p'{p)/p > 0, for p > 0, (3.46) 

so that 

j A{Tp')dx > j A{p)dx + j A!{p){Tp' - p), for i > L (3.47) 
bmce A'{p) E and Tp* weakly converges to p in L'^ , 

J A'{p)(Tp' - p) ^ 0, as i ^ +00. (3.48) 



Therefore, 



[ A{p)dx < lim inf / A{Tp')dx. (3.49) 
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Next, following the proof in [28], we can show that 



, Tp^(x)L(mrp^i(r(x)) — p(x)L(mo(r(x)) 
lim inf / ' ^ ^ ' ' ^ ^ ^' dx > 0, 3.50 



t— >oo 



by showing that the mass function 

p{x)dx 



is continuous for r > 0, and using (j3.6|) . Then (j3.44p follows from (j3.43p . (j3.49p and (|3.5Up . 
3.2 Stability 

In this section, we assume that the pressure function p satisfies 

peC^[0,+oo), lim47F = 0, lim ^ = K, > for p > 0. (3.51) 

where < K < +oo and 7 > 4/3 are constants. It should be noticed that (|3.5ip implies both 
(2.15) and (|3.3p . We consider the Cauchy problem for (jl.ip with the initial data 

p{x, 0) = po{x), ^r{x, 0) = vo(x). (3.52) 

We begin by giving the definition of a weak solution. 

Definition: Let pv = m. The triple {p,m,^){x,t) (x £R^,t e [0,T]) (T > 0) and <I> given 
by (|1.2p . with p > 0, m, ui0ui/p and pV^ being in Lj^^{M.'^ x [0, T]), is called a weak solution 
of the Cauchy problem (jl.ip and (j3.52p on x [0, T] if for any Lipschitz continuous test 
functions ip and ^ = (^/'i, "02; V^s) with compact supports in x [0, T], 

T 

' (pV't + m • Vip) dxdt + / pq{x)iIj{x, Q)dx = 0, (3.53) 





and 

j j • + • V^^ dxdt + y mo(a;)*(x,0)(ix = ^ ^ pV^^dxdt, (3.54) 

both hold. 

For any weak solution, it is easy to verify that the total mass is conserved by using a 
generalized divergence theorem for L'^ functions (r > 1) (cf. p]), 

j p{x,t)dx = j p{x,0)dx, t>0. (3.55) 

The total energy of system (jl.ip at time t is 

E{t) = E{p{t)Mt)) = j I^Mp) + ^P|v|') {x,t)dx -^J \V<^\\x,t)dx, (3.56) 
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where as before, 

A{p) = P r ^ds. (3.57) 
Jo 

For a solution of (|l.ip without shock waves, the total energy is conserved, i.e., E(t) = E{0) 
{t > 0)(cf. [35j). For solutions with shock waves, the energy should be non- increasing in 
time, so that for all t > 0, 

E{t) < ^(0), (3.58) 

due to the entropy conditions, which are motivated by the second law of thermodynamics (cf. 
[20] and [34]). This was proved in [28]. 

We consider axi-symmetric initial data, which takes the form 

Po{x) = p(.r,z), 

vo{x) = f5(r, z)er + Uo(r, z)eg + v^{p, z)e3. (3.59) 
Here r = ^yxf + X2, z = X3, x = {xi,X2,X3) G (as before), and 

e, = (xi/r, X2/r, of, eg = (-xa/r, xi/r, of, eg = (0, 0, 1)^. (3.60) 
We seek axi-symmetric solutions of the form 

p{x,t) = p{r,z,t), 

v(x, t) = v''{r, z, t)er + v^{r, z, t)eg + v^{r, z, t)e3, (3.61) 
q>[x, t) = ^>(r, z, t) = -Bp{r, z, t), (3.62) 

We call a vector field u{x,t) = {ui,U2,U3){x) (x G M'^ ) axi-symmetric if it can be written in 
the form 

u(x) = u^{r, z)er + u^{r, z)eg + u^{p, z)e3. 

For the velocity field v = {vi,V2,V3){x,t), we define the angular momentum (per unit mass) 
j{x,t) about the X3-axis at {x,t) , t > 0, by 

j{x, t) = X1V2 - X2V1. (3.63) 

For an axi-symmetric velocity field 

v(x, t) = v'^{r, z, t)er + v^{r, z, t)eg + v^{p, z, t)e3, (3.64) 
vi = — V V , V2 = — V H V ,V3 = V , (3.65) 

so that 

j{x,t) =rv\r,z,t). (3.66) 
In view of ( I3.64p and (j3.66p , we have 

|v|2 = It;^^ ^ ^ + |t>3|2. (3.67) 
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Therefore, the total energy at time t can be written as 

\yBpWx,t)dx + ^ J p{\v^\'' + \v^\''){x,t)dx. (3.68) 

There are two important conserved quantities for the Euler-Poisson equations (jl.ip : namely 
the total mass and the angular momentum. In order to describe these, we define Dt, the non- 
vacuum region at time i > of the solution by 

Dt = {x : p{x, t) > 0}. (3.69) 

We will make the following physically reasonable assumptions Al)-A4) on weak solutions of 
the Cauchy problem and (|3.52|) (Al)-A4) are easily verified for smooth solutions. For 
general weak solutions, they are motivated by physical considerations, cf. |35j ) . 

Al) For any t > 0, there exists a measurable subset Gt C Dt with meas{Dt — Gt) = 
(meas denotes Lebsegue measure) such that, for any x £ Gt, there exists a unique (backwards) 
particle path ^{T,x,t) for < r < t satisfying 

dr^{T, X, t) = v(e(T, X, t), r), i{t, X, t) = X. (3.70) 

Remark 3. If v(-,t) E BV{^^) and div^v{-,t) G L°^(M3) for t > {div^ is in the sense of 
distributions), it was proved by L. Ambrosio ([!]) that Al) is valid. Related results can be 
found in [9]. 

For X £ Gt, we write 

c(o,x,t) = e-t(x). 

Also, for X G and t > 0, we denote the total mass at time t in the cylinder {y G : 
r{y) < r{x)} by m^(t)(r(x)), i.e., 

mp{t)ir{x)) = p{y,t)dy. (3.71) 

Jr{y)<r{x) 

For axi-symmetric motion, we assume 
A2) 

mpit){r{x)) = mpo(r(e-t(x))), for xeGt,t>0. (3.72) 

(This means that the mass enclosed within any material volume cannot change as we follow 
the volume in its motion ( [35j, p. 47)). Moreover, we assume that the angular momentum 
is conserved along the particle path: 

A3) 

jix,t) = jiC-tix),0), iovxeGt,t>0. (3.73) 
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Finally, for L = j^, we need a technical assumption; namely, 
A4) 

lim ^ "'^''^ ' " "^''^ ' = 0, (3.74) 

for i > 0, where (j{t) = p{t) — p. 

Remark 4. (j3.74p can be understood as follows. For any p G Wm, we have limr^o+ iT^pif) = 0. 
Therefore limr^o+ -^("ip(t)(^) + f^p{^)) = ^(0) = 0, so if we define 

/+00 
(/)(s,z,t) - p{s,z))dz, 
-oo 

then if 

^<r(t){r) J^{2TTs{p{s,t)- ^{s))ds 

^5 = — K G L [0,0) for some > 0, (3.75) 

(|3.74p will hold, li p{- , t) —~p{-) G L°°(0, 5), then (j3.75p holds. This can be assured by assuming 
that p{r, z, t) — p{r, z) G L°°((0, (5) x M x M+) and decays fast enough in the z direction. For 
example, when p{x,t) — p{x) has compact support in and p{-,t) — p(-) G L°°(M'^), then 
([3:71]) holds. 

We next make some assumptions on the initial data; namely, we assume that the initial 
data is such that the initial total mass and angular momentum are the same as those of the 
rotating star solution (those two quantities are conserved quantities). Therefore, we require 

Ii) 

j po{x)dx = j p{x)dx = M. (3.76) 

Moreover we assume 

I2) For the initial angular momentum j(x, 0) = rvQ{r, z) =: jo{r, z) (r = y^a;^ + x^, z = x^, 
for x = (xi,X2,X3), we assume j(^,0) only depends on the total mass in the cylinder {y G 
'M?,r{y) < r{x)}, i.e. , 

j{x,0)=jo{mpMx)))- (3.77) 
(This implies that we require that VQ^r, z) only depends on r.) 

Finally, we assume that the initial profile of the angular momentum per unit mass is the same 
as that of the rotating star solution, i. e., 

I3) 

jl{m)=L{m), 0<m<M, (3.78) 

where L{m) is the profile of the square of the angular momentum of the rotating star defined 
in Section 2. 
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In order to state our stability result, we need some notation. Let A be the constant in Theorem 
2.2, i.e., 

+ /~ ) L{mpis)).s-^ds - Bp{x) = A, x G T, 
l^^^L{mp){s))s~^ds-Bp{x)>X, xeR^-F, 

with A defined in (13.57P and F defined in (2.16). 

For p £ Wm, we define, 

dip, p)= [ [Aip) - Aip)] + ip-p) r {^^^^^^ds - A - Bp}dx. (3.80) 

J Jr{x) S'^ 

For X G r, in view of the convexity of the function A (cf. (j3.46p ) and (j3.79p . we have, 
(Aip) - A{p)){x) + ( r ^^"^i'^^ ds - A - Bp{x)){p - p) 



J r{x) ^ 

= {A{p) - A{p) - A'{p){p - p)){x) > 0. (3.81) 

For X G — r, p{x) = 0, so we have A{p){x)) = 0. This is because since ^(0) = due to 
p{0) = (cf. (3.3)) and (2.5). Therefore, by UjSTld^i . we have, for p G Wm and x G - F, 

{Aip) - ^(/5))(x) + ( r ^^"^f^^ s - A - Bpix))ip - p) 

J r{x) ^ 

= Aip) > 0. (3.82) 
Thus, for p G Wm, 

dip,p)>0. (3.83) 



We also define 



, / ~N 1 /■ p(x)L(mp(r(x)) - p(x)L(mp(r(x)) 
diiP^P) = 2 J ;Ji^) 



oo 

s-^Limpis))dsipix) - pix))dx, (3.84) 

r{x) 

for p G Wm- We shall show later that di > 0. Our main stability result in this paper is the 
following global-in-time stability theorem. 



Theorem 3.2. Suppose that the pressure function satisfies 13. 51\) . and both ^3.4\ ), (37l 
hold. Let p be a minimizer of the functional F in Wm, and assume that it is unique up to a 
vertical shift. Assume that Ii)- I^), 70^ - \3. 78^ ] hold. Moreover, assume that the angular 
momentum of the rotating star solution p satisfies 12. 9\) . h3. 1\) and (3.2). Let (p, v, <l>)(x, t) 
be an axi-symmetric weak solution of the Cauchy problem and \3. 52\) satisfying the 
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assumptions Al)-A4), [ \3. 70[ j-( fg. 74\) ]- If the total energy E{t) (of. 113. 56\) ) is non-increasing 
with respect to t, then for every e > 0, there exists a number (5 > such that if 



d{po,p) + ^WVBpo - VBpWl + \di{po,p)\ 



+ 



/Oo(x)(|fSP + \v^\'^){x)dx < S, 



3|2n 



(3.85) 



then there is a vertical shift ae^ (a G M, es = (0, 0, 1) ) such that, for every t > 
d{p{t),T-p) + j^W^Bpit) - VBT-p\\l + |di(p(t),r-^)| 



+ 



p{x,t){\v''{x,t)\'^ + \v^{x,t)\'^)dx < e, 



(3.86) 



where T"'p{x) =: p{x + aes). 



Remark 5. As noted in |28J, the vertical shift aes appearing in the theorem is analogous to 
a similar phenomenon which appears in the study of stability of viscous traveling waves in 
conservation laws, whereby convergence is to a "shift" of the original traveling wave. 

Remark 6. Without the uniqueness assumption for the minimizer of F in Wm^ we can have 
the following type of stability result, as observed in ^2] for the non-rotating star solutions. 
Suppose the assumptions in Theorem 13.21 hold. Let Sm be the set of all minimizers of F in 
Wm and (p, v, $)(x,f) be an axi-symmetric weak solution of the Cauchy problem (jl.ip and 
()3.52p . If the total energy E{t) is non-increasing with respect to t, then for every e > 0, there 
exists a number 5 > such that if 



inf 

P&Sm 



d{po,p) + ^llVBpo - ^Bp\\l + \di{po,p)\ 
ovr 



+ 2 / Po{x){Kr + \v^\'){x)dx <5, 



(3.87) 



then for every t > 



inf 

P&Sm 

+ ■ 



d{p{t),T-p) + ^\\^Bp{t) - VBT'^pWl + \d^{p{t),T-p)\ 
p{x,t){\v''{x,t)\^ + \v^{x,t)\^){x)dx < e. 



then there is a vertical shift aes (a G M, es = (0, 0, 1)) such that, for every t > 
d{p{t),T-p) + ^l|Vi?p(0 - VBT-pWl + \d^{p{t),T'^p)\ 



+ 



p(x,t)(|t;'^(x,t)p + \v^{x,t)\'^)dx < e. 



(3. 



(3.89) 



where T^p{x) =: p{x + aes). 
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The proof of Theorem 13.21 follows from several lemmas. The proofs of these lemmas are 
similar to those in [28], and therefore we only sketch them. First we have 

Lemma 3.8. Suppose the angular momentum of the rotating star solutions satisfies i2. 9\) . 
/ TO]] and (3.2). For any p{x) € Wm, if 

lim L{mp{r) + mp{r))mcj[r)r~'^ = {), (3.90) 

where a = p — p, then 

di{p,p)>0, (3.91) 

where di is defined by 1^3.84^ . 

Proof. For an axi-symmetric function f{x) = f{r,z) (r = y^xf + z = X3 for x = 
(xi,X2,X3)), we let 

r+00 

/(r)=27rr/ f{r,z)dz, (3.92) 



m/(r)= / f{x)dx= / f{s)ds, (3.93) 

J{x:^/xl+x^<r} Jo 

SO that 

m'f{r) = f{r). (3.94) 

In order to show (j3.9ip . we let 

a{x) = {p-p){x), (3.95) 

and for < a < 1, we define 

If {p + aa){x)L{mp+a^{r{x))) - p{x)L{mp{r{x))) 

Q(-) = 2 J ^) 

p poo 

— a / s^^ L{mp{s))dsa{x)dx. (3.96) 

J Jr{x) 

Then 

Q(0) = 0, Q(l) =di(/>, ~p). (3.97) 
By the assumption that V{m) > for < m < M (cf. (3.2)) and (|3.90p . we can show that 

Q\a) = a{r) J s-^{L{mp+^^{s)) - L{mp{s)))dsdr, (3.98) 

and therefore 

Q(0) = Q'(0) = 0. (3.99) 

This is done by interchanging the order of integration and integrating by parts (details can 
be found in [28j). Differentiating (j3.99p again and interchanging the order of integration, we 
get 

d^Q{a) 



-«/ s-'L'{mp+^^{s)){m„{s)yds. (3.100) 



21 



Therefore, if L'{m) > for < m < M, then 

> 0, /or < a < 1. (3.101) 
This, together with ([3:99]) and ([SJZ]), yields di{p,p) = Q(l) > 0. □ 

Lemma 3.9. Let (p, v) be a solution of the Cauchy problem U.l\) . 13. 52\) as stated in Theorem 
3.2, then 

E{p,^){t)-F{~p) 

= dipit),p)+diipit),p) - ^\\VBpi;t) - VBpWl 

+ \j + (3.102) 

Proof. From Al)-A3), we can show 

j2(x,t) =L(m^(t)(r(x))), x d Gt. (3.103) 

Therefore, by (j3.68|) . we have 

/I f p(x,t)L(m„(A(r(x)) 
A{p)ix,t)dx + -J ^-^^JA^^^^^Jidx 

j \VBp\^{x, t)dx + pi\v'\'^ + \v^\'^)ix, t)dx. (3.104) 
(|3.102l) fohows from (I3.104j) and the follow identities: 
{\\VBp{.Ml-\WBp\\l) 

= \ \V{Bp{-,t)) - VBp)\\l + 2 j VBp{x) ■ {VBp{x,t) - VBp{x))dx 
= \\V{Bp{;t))-VBp)\\l-87r J Bp{x){p{x,t) - p{x))dx. 

and 



fpix,t)dx = jp{x)dx = M. 



□ 



Having established these lemmas, the proof of Theorem 3.2 is the same as the proof of 
Theorem 3.1 in 128). 
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4 Applications to White Dwarf and Supermassive Stars 



In this section, we want to verify the assumptions (j3.4p and (jS.Sp in Theorem 3.1 for both 
white dwarfs and supermassive stars. Once we verify (j3.4p and (j3.5p . we can apply Theorems 
3.1 and 3.2. We begin with the following theorem which verifies (3.5) for white dwarfs, 
supermassive stars, and polytropes with 7 > 4/3, in both the rotating and non-rotating 
cases. 

Theorem 4.1. Assume that the pressure function p satisfies (3.3). Then there exists a 
constant Tic satisfying < Tic < 00 if ^ = 4/3 and Tic = 00 if ^ > 4/3, such that if 
M <Tlc, then [3~5\) holds for p £ Wm- 

Proof. Using (j3.13p . we have, for p G Wm, 

J 2 r{xY J 

Taking p = 1, (7 = 4/3, r = 7, and a = (where 7 > 4/3 is the constant in (3.3)) in 

Young's inequality (j3.1ip . we obtain, 

IIP||4/3<IIP||?IHI^-" = M'^IMI^-"- (4.2) 
This, together with (3.16)-(3.18) yields 

j p^'^dx < Mi"( j p^dxf < mI" (^{p*y-^M + aj A{p)dx^ 

< C (^Mt"+''(/9*)^/3 ^ aM^i j A{p)dx)^^ , (4.3) 

where b = g^^^-^^ , a and p* are the constants in (3.17) and we have used the elementary 
inequality (x + y)^ < C{x^ + y^), for x, y > 0, < 6 < 1, for some constant C . Therefore, 
(f4T]l and (f43]l imply 

j + < F{p)+CaMh+l( [ A{p)dxf+CM^'^+'+l{p*f'\ (4.4) 

J 2 r{xY J 

If 7 > 4/3, then < 6 < 1, if 7 = 4/3, then 6 = 1. Therefore dO]) implies ([33]). □ 

The next result shows that (3.4) holds for a wide class of (rotating or non-rotating) stars, 
including White Dwarfs. 
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Theorem 4.2. Suppose that the pressure function p satisfies 13. 3\) and 

lim ^ = /?, (4.5) 

for some constants 71 > 4/3 anc? < /3 < +cxo, and assume that the angular momentum (per 
unit mass) satisfies (2.9). Then there exists satisfying < < +00 if ^ = 4/3 and 
Mc = +00 i/7 > 4/3 such that if M < M.^, then holds, where 7 is the constant in i3.3\) . 

Remark 7. White dwarfs satisfy (13. 3p and (14. Sp with 7 = 4/3 and 71 = 5/3. 
Proof of Theorem \4-2\ 

Due to (j3.3p and ()4.5p . we can apply Theorem 2.1. Let /5(x) G Wm,5 be a minimizer F(/5) in 
Wm,s as described in Theorem 2.1, and let 

G = {xGR^ : p{x) > 0}. 

Then G is a compact set in M^, and p € C^{G). Furthermore, there exists a constant ^ < 
such that 

A'ipix)) + I^^)Hmpis)s-^ds-Bp{x)=p, x€G, ^^^^ 

f^^^ L{mf,{s)s~Hs - Bp{x) >is, xeR^-G. 

It follows from [2j that there exists p G Wm,s C Wm such that F{p) = infpgpi/j^,^ g F{p)- It 
is easy to verify that the triple (p, v, <1>) is a time-independent solution of the Euler-Poisson 
equations <^ in the region G = {x G . > g}, where v = (_ ^2J(mp(r)) ^ xiJ(m^(r)) ^ 

and $ = —Bp. Therefore 

V:,p(/5) = pV^{Bp) + pL{mfi)r{x)-^er, x G G, (4.7) 

where e,. = (;^, 0). Moreover, it is proved in [1] that the boundary dG of G is smooth 
enough to apply the Gauss-Green formula on G. Applying the Gauss-Green formula on G 
and noting that p\qg = 0, we obtain. 



X ■ Vxp{p)dx = —3 / p{p)dx = —3 / p{p)dx. (4-8) 
G Jg J 

As in [28], we have 

/ X ■ pV xBpdx = — / pBpdx = — j pBpdx. (4-9) 
JG Jg 27 

Next, since x ■ er = r{x), we have 

X ■ p{x)L{mp{r{x))r^^{x)erdx 

G 

p{x)L{mp{r{x))r^'^ {x)dx 

G 

p{x)L{mp{r{x))r^'^{x)dx. (4-10) 
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Therefore, from (fi^ - ()4.10p we have 

~ J Pip)dx = J pBpdx + j p{x)L{mp{r{x))r^'^{x)dx. (4-11) 

Let p{x) = b^p{bx), for 6 > 0; then p G Wm- Also, it is easy to verify that the following 
identities hold, 

-pB-pdx = f [ P\^dxdy 



x-y\ 

x-y\ 

A{p)dx = b-^ I A{b-'p{x))dx. (4.13) 



b I I I P\^}P^dxdy = b I pBpdx (4.12) 



Moreover, for r > 0, 



r 



m p{r) = 27r s p{s,z)dsdz 
Jo J-oo 
f^r poo 

= 2tt I si p{bs,bz)dsdz 

Jo J-oo 
i-br foo 

= 2tt s' p{s',z')ds'dz' 

Jo J~oo 



Therefore, 



mp(br). (4.14) 



p{x)L{mp{r{x))) f b^ p{x)L{m^{br{x))) 



dx = I r-rr; dx 

dx. (4.15) 



r(x)^ J r{x) 

_,2 f p{x)L{mi,{r{x))) 



r{xy 



It follows from (flJ2]) - (jlJ5]) that 

F{p) = j A{b^p)dx j pBpdx 



+ ^ / P(^)^K(^(^)))rf,. (4.16) 
2 J r{xY 



Hence, diH]) and (|4l6l) give 

F{p) = [ {b-^A{b^p) - Sbpipix))) dx 



^--b) [ Pm-ip-Rdx. (4.17) 
2 J J r[xY 



In view of (2.9), we have 

'6^ _ \ r p{x)L{mp {r{x))) 
2 ) j r{x) 



"--b] j ' "' dx < 0, (4.18) 
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if & > is small. It follows from (13.91) that 



^Pp'" < pip) < 2pp'^\ for small p. (4.19) 



Thus, when b is small, since p is bounded, we have 



^ 6371 (^)7i(^) < A{b^p{x)) < ^^6371(^)71(3;)^ (4.2O) 



2(71 - 1) vr, V . - V rv - 

for X e M^. Hence, KT8\\ and (f4T9]l imply 

[b-'^^Aib^p) - 3bp(p{x))) dx 



Since 71 > 4/3, we have 371 — 3 > 1. Therefore, we conclude that 

[b-^A{b^p) - 3bp{p{x))) dx < 0, (4.22) 
for small b. (j3.4p follows from (j4.17p , (j4.18p and ()4.22p . This completes the proof of Theorem 

4.2.n 

We show next that if the angular momentum distribution is everywhere positive, we may 
apply the existence theorem of Friedman and Tarkington, [13], to conclude that (3.4) holds 
with no total mass restriction. This result applies also to White Dwarfs. 

Theorem 4.3. Suppose that the pressure function p satisfies i3. 3\) with 7 = 4/3 and \3. 9\) 
holds. Assume that the angular momentum (per unit mass) J{m) = L{m) satisfies (2.14), 
then |g.^[) holds for p G Wm for any < M < +00. 



Proof. By the existence theorem in [13j . if (2.14) is satisfied, then for any < M < +00, there 
exits p G Wm,s such that F{p) = inip^wM s ^(p)- Also, all the properties of p in Theorem 
2.1 are satisfied. Moreover, the regularity of the boundary dG is smooth enough to apply the 
Gauss-Green formula (cf. [4]). The proof now follows exactly as in Theorem 4.2. □ 

We finally turn to the case of rotating supermassive stars. 

Theorem 4.4. Consider suppermasive star; i.e., 

p{p) = kp^l'^ , A; > is a constant. (4.23) 



If there exists p G Wm such that p G C^{G) n C(M^) is a steady state solution of the Euler- 
Poisson equation with the velocity field v = (_ ^2v^(m-p(r)) ^ xi%/L(mp{j)) ^ open bounded 

set G C with the Lipschitz boundary dG, i.e. 



^xP{p) = P^xiBp) + pL{mp)r{x) ^e^, x G G, 

(4.24j 

p = 0, x G - G. 
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then holds provided L satisfies (2.9) and 

L{mo) > 0, for some mo G (0,M). (4.25) 

Proof. Following along the same lines as (4.7)-(4.10), we obtain the same equality as (4.11). 
Therefore, 

F{p) = -^J p{x)L{mp{r{x))r-^{x)dx, (4.26) 

in view of and (4.11). Since p e C^{G) D C{R^) and /5 = for x e - G, it is 

easy to show that mp{r) is continuous in r. Moreover, m^(0) = and mp{R) = M, where 
R — max2,gQ(r(x) . Therefore, there exits tq G (0,M) such that 

"mpiro) = mo, (4.27) 

where mo is the constant in (I4.25p . Thus, 

L{mp{ro)) > 0, (4.28) 

in view of ()4.25p . Since m^{r) is continuous in r and L{m) is continuous in m, we conclude 
that 

p{x)L{mp{r{x))r~'^ {x)dx > 0. (4.29) 

The inequality (j3.4p now follows from (j4.26p ). □ 

The preceding theorems, together with Theorem 3.2 show that polytropes {p{p) = kp^) 
with 7 > 4/3 and White Dwarf stars, in both the rotating and non-rotating cases, as well 
as rotating Supermassive stars are dynamically stable. Moreover, if the angular momentum 
distribution is not everywhere positive and the pressure p behaves asymptotically near infinity 
like then dynamic stability holds only under a (Chandrasekhar) mass restriction, M < 
Mc. 

5 Nonlinear Dynamical Stability of Non-Rotating White Dwarf 
Stars With General Perturbations 

The dynamical stability results in Section 3 apply for axi-symmetric perturbations. In this 
section, we prove the nonlinear dynamical stability for non-rotating white dwarf stars with 
general perturbations. For white dwarf stars, as mentioned before, the pressure function 
satisfies 

p e C^[0,+cx)), lim^ = /3, lim ^ = i^:, p'(p) > for p > 0, (5.1) 

p^0+ p'l'i p^oo p"/ 
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where 71 > 4/3, < /3 < +00 and < K < +00 are constants. In this section, we always 
assume that the pressure function satisfies (j5.ip . First, we define for < M < +00, 



Xm = {/o : ^ M,p > 0,o.e., p{x)dx = M. 



I 



[A{p{x)) + ^p{x)Bp{x)]dx < +00}, (5.2) 



where A{p) is the function given in (2.5). For p £ Xm, we define the energy functional G 
for non-rotating stars by 

G{p) = I [A{p{x)) - ^pix)Bpix)]dx. (5.3) 
We begin with the following theorem. 

Theorem 5.1. Suppose that the pressure function p satisfies i5.1\) . Let p^ be a minimizer 
of the energy functional G in Xm and let 

Tn = {x£R^ : pn{x) > 0}, (5.4) 

then there exists a constant Xn such that 

A'{pn{x)) - Bpm{x) = \n, x e Fat, 

(5.5) 

-Bpn{x)>\n, xSM^-Fat. 

The proof of this theorem is well-known, cf. [32] or [2]. 

Remark 8. 1) We call the minimizer pn oi the functional G in Xm a non-rotating star 
solution. 

2) It follows from [23] that the minimizer pN of the functional G in Xm is actually radial, 
and has a compact support. 

Similar to Theorem 3.1, we have the following compactness theorem. 

Theorem 5.2. Suppose that the pressure function p satisfies i5.1\) . There exists a constant 
M" (0 < JVr < 00) such that if M < M", then the following hold: 

(1) 

inf G{p) < 0, (5.6) 

P&Xm 

(2) for p e Xm, 

j A{p){x)dx < CiG{p) + G2, (5.7) 
for some positive constants Ci and G2, 

(3) if{p^} C Xm is a minimizing sequence for the functional G, then there exist a sequence of 
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translations {x"^} C M^, a subsequence of {p^}, (still labeled {p^}), and a function pN G Xm, 
such that for any e > there exists R > with 

Tp\x)dx <e, ien, (5.8) 

\x\>R 

and 

Tp\x) pN, weakly in L^/^(M^), as i ^ oo, (5.9) 

where Tp^{x) := p^{x + x^). 
Moreover 

(4) 

VB{Tp') VB{pn) strongly in L^{W^), as i ^ oo, (5.10) 

and 

(5) p is a minimizer of G in Xm- 

Proof. First, the proofs of (1) and (2) are the same as Theorems 4.1 and 4.2 by taking L = 
(it is easy to check the axial symmetry is not used the the proof of Theorems 4.1 and 4.2 
if L = 0). Lemmas 3.4, 3.5 and 3.7 still hold by taking 7 = 4/3 and L = 0, and replacing 
Wm by Xm, F by G and /m by infpg^M G{p). Also, it is easy to check that (3.25)-(3.29) in 
the proof of Lemma 3.6 still hold by replacing fM by inf pg^M G{p). Therefore, following the 
proof of Lemma 3.6, we conclude: 

If C Xm is a minimizing sequence for G, then there exists constant 5o > 0, io G N and 
X* G M^, such that 

p^{x)dx > (5o, i > io- 



Therefore, if we let 



then 



Tp\x) := p\x + x'), (5.11) 



/ Tp'^{x)dx > 5q, i > if). 
Jbi{0) 



This is similar to (3.39). Having established this inequality and the other analogues of 
Lemmas 3.4, 3.5 and 3.7, we can prove this theorem in a similar manner as the proof of 
Theorem 3.1. □ 

For the stability, we consider the Cauchy problem (1.1) with the initial data (3.53). We 
do not assume that the initial data have any symmetry. 

Let pn he a minimizer of G on Xm and Aat be the constant in (|5.5p . For p £ Xm, we define 



d{p, pn) = / {[A{p) - A{pn)] - {p- Pn){>>n + BpN}dx, 

{[Aip) - A{pn)] - BpNip - PN)]dx, (5.12) 
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where we have used the identity 

J pdx = J p^dx = M, 
for p £ Xm- By a similar argument as (3.82) and (3.83), we have 

d{p,pN)>Q, (5.13) 

for any p G Xm, in view of (j4.6p . Our nonhnear stabihty theorem of non-rotating white dwarf 
star solutions is the following theorem, which extends the results in |32j. 

Theorem 5.3. Suppose that the pressure function satisfies h5.1\) . Let p^ he a minimizer of 
the functional G in Xm, (ind assume that it is unique up to a translation pn{x) pj\f{x + y). 
Let (/9, V, <I>)(x, t) be a weak solution of the Cauchy problem II 1-1]) and liS. 52\) satisfying 



p{x,t) = J po{x) = J pN{x)dx = M. 

If the total energy E{t) ( cf. h3. 56\) ) is non-increasing with respect to t, then there exists a 
constant M'^ (0 < M'^ < oo) such that if M < M^, then for every e > 0, there exists a number 
6 > such that if 

d{po,PN) + -^W^Bpo - VBpnWI + po{x)i\vo\'^)ix)dx < 5, (5.14) 
then there is a translation y G such that, for every t > 

d{p{t),TypN) + ^\\^Bp{t) - VBTy-pNWl + \j p{x,t)\v{x,t)\'')dx < e, (5.15) 

where T^pN^x) =: pn{x + y). 

The proof of this theorem follows from the compactness result (Theorem 15. 2p , and the 
arguments in [28] and [32], and is thus omitted. 
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